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Cryogenic propellants LOx/H2 are used at very high pressure in rocket engine combustion. The description of the 
combustion process in such application is very complex due essentially to the supercritical regime. Ideal gas law 
becomes invalid. In order to try to capture the average characteristics of this combustion process, numerical 
computations are performed using a model based on a one-phase multi-component approach. Such work requires 
fluid properties and a correct definition of the mixture behavior generally described by cubic equations of state 
with appropriated thermodynamic relations validated against the NIST data. In this study we consider an alterna- 
tive way to get the effect of real gas by testing the volume-weighted-mixing-law with association of the compo- 
nent transport properties using directly the NIST library data fitting including the supercritical regime range. The 
numerical simulations are carried out using 3D RANS approach associated with two tested turbulence models, the 
standard k-Epsilon model and the realizable k-Epsilon one. The combustion model is also associated with two 
chemical reaction mechanisms. The first one is a one-step generic chemical reaction and the second one is a 


two-step chemical reaction. The obtained results like temperature profiles, recirculation zones, visible flame 


lengths and distributions of OH species are discussed. 
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Introduction 


Rocket-engine based on cryogenic propellant tech- 
nology has been used for many years in launchers like 
those of Ariane 1 to 5. The use of liquid propellant allows 
significant progress due to its high performance and good 
flexibility. Cryogenic hydrogen as fuel and liquid oxygen 
as oxidizer are injected separately into the combustion 
chamber at a very high pressure beyond the thermody- 
namic critical point. This mixture condition is well 
known as an interesting solution that provides best effi- 
ciency and performance for a variety of launch applica- 
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tions [1]. However, reaching equilibrium between low- 
cost and short-time development process remains a chal- 
lenging issue. The degree of local mixing and the extent 
of the combustion reactions must be established in order 
to accurately predict the performances and the internal 
thermal environment of the engine. The lack of experi- 
ments to reproduce in the laboratory environment the 
conditions existing in the combustion chamber and the 
limited scientific knowledge in this field lead to numeri- 
cal studies which are currently performed by space agen- 
cies. Some of these numerical simulations based on LES, 
DNS, or RANS approaches were detailed in [2]. For 
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Nomenclature 

Cp Cy Specific heat capacity at constant pressure 
and volume, J.kg.K 

D Diameter of the injector, m 

D, Diffusion coefficient of a component a, 
m’.s™ 

e Internal energy, J.kg ` 

ha Enthalpy of the component, J.kg ' 

k Turbulent kinetic energy, m’.s* 

I Identity tensor 

M Molecular weight, g.mol' 

Pr Prandtl number 

p Static pressure, Pa 

q Heat flux, W.m” 

Se Schmidt number 

T Temperature, K 

Ui Velocity, i" component, m.s" 

V Molar volume, m*.mol”! 

Y, Mass fraction of a component a 

Greek symbols 

E Dissipation rate of k, m°.s” 

À Thermal conductivity, W.m K 

u Dynamic viscosity, kg.m'.s 

Oy Net production rate of a component, 
kg.m?.s | 

p Density, kg.m™ 

T Viscous stress tensor 

Subcripts 

c Critical 


example, five different simulations, grouped in the pro- 
ject supported by the NASA agency and called Constel- 
lation [3] were compared to extract a global behavior. 
The results of this study point that the use of the LES or 
RANS codes shows no similarity in the results. One rea- 
son is due to subgrid models used for high pressure ap- 
plications as they were originally developed for cases at 
low pressure. Other studies have been conducted at the 
European research center CERFACS by Thomas Schmitt 
[4]. The first tests therefore show the importance of bet- 
ter understanding the various unknowns of the problem. 
Even if LOx-H, rocket engines have been relatively safe 
during the last years, the processes that control combus- 
tion are still not well understood [5]. For instance, both 
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F Formation 

t Turbulent 

a Reynolds averaging 

a Favre averaging 

j Fluctuation 

Notations 

AMR Adaptive Mesh Refinement 

BWR Benedict-Webb-Rubin 

CERFACS European Centre for Research and Ad- 
vanced Training in Scientific Computation 

CFD Computational Fluid Dynamics 

DNS Direct Numerical Simulation 

EDM Eddy Dissipation Model 

GH2 Gaseous Hydrogen 

IWRCM International Workshop on Rocket Com- 
bustion Modelling 

Ke Standard (k-e) turbulence model 

LES Large Eddy Simulation 

LOX Liquid Oxygen 

LRE Liquid Rocket Engine 

NIST National Institute of Standards and Tech- 
nology 

ONERA Office National d’Etudes et de Recherches 
Aérospatiales 

PR Peng-Robinson state equation 

RKE Realizable (k-€) turbulence model 

SRK Soave-Redlich Kwong equation 

VWML Volume Weigthed Mixing Law 


for surrounding quiescent or moving gases and for su- 
percritical shear flows such as those occurring in the 
rocket combustors, research works are in progress [6-9]. 
For one decade, a series of investigations have been per- 
formed based on a research-test facilities with high pres- 
sure combustion of cryogenic propellants, namely with 
the well known Mascotte device (ONERA, France) and 
the European research and technology test bench P8 
[10,11,12-17]. With the advantage of optical access to the 
combustion chamber, experimental visualization using 
advanced laser techniques are possible mainly like Co- 
herent Anti-stokes Raman Spectroscopy (CARS), Planar 
Laser Induced Fluorescence (PLIF) and line-of-sight 1m- 
aging (shadowgraphs and OH* emission). Also, accurate 


Abdelkrim Benmansour et al. Study of LO2/GH2 Supercritical Combustion in the ONERA- Mascotte Test-rig Configuration 99 


monitoring of LOx/GH, flame shapes and lengths is al- 
lowed, as well as an identification of controlling factors. 
In the same way, several numerical investigations were 
done for turbulent mixing with crossing the critical fluid 
pressure or temperature [18,19]. Schmitt et al. investi- 
gated the effect of chamber pressure on the turbulent 
mixing LOx /GH, and examined the effect of injection 
temperature [20]. The stabilization of the LOx/GH) flame 
tip at the injector rim was investigated in [21-22]. Lami- 
nar counter flow diffusion flames have been studied to 
investigate the chemical structure of LOx/GH), transcriti- 
cal flame [23-26]. The coaxial jet flame was investigated 
using a RANS approach in [27-29]. Other researchers 
have performed numerical investigations to try to depict 
the complex mixing and combustion processes occurring 
in rocket combustion chambers. The main problem is to 
get the behavior of the propellant mixture at the super- 
critical conditions. Indeed, the evaluation of thermody- 
namic properties and transport parameters constitute the 
major difficulty to describe correctly the mixing and 
combustion processes. Since the year 2005, some studies 
[30-33] have used the LES approach to improve the un- 
derstanding of fluid dynamics inside LREs. So, modern 
theory, however, still lacks a physically-based model to 
explain the phenomenology when fluid exceeds the 
thermodynamical critical condition [34]. For instance, in 
references [35-41] studies are made on the behavior, 
modeling and analysis of turbulent supercritical mixing. 
Reference [40] performed numerical investigations using 
DNS approach to study the temporal non-reacting mixing 
layers, the impact of real-gas thermodynamics on the 
characteristics of turbulent mixing and to determine the 
validity of closure terms for Large Eddy Simulation 
(LES). However, nowadays these methods are still too 
computationally expensive to be used in industrial appli- 
cations knowing that it is currently impossible to achieve 
simulations identical to those of the Constellation project 
[3]. Currently, only DNS to our knowledge, which is ap- 
plied to cryogenic engines, was carried out by Bellan’s 
team [40] and F. Troadec [2] where the density ratio of 
the mixture (O2/H2) hardly exceeds 26 and which is still 
very far from real applications where it is greater than 
100. 

This work deals with a numerical study carried out on 
the combustion processes occurring in the Mascotte test 
rig (V03). A 3D RANS study of LO /GH) supercritical 
combustion is accomplished. Real gas effects based on 
polynomial functions for transport parameters and densi- 
ty of pure species have been introduced into the com- 
mercial code Fluent through user defined functions. The 
homogenous multi-component, one-phase model asso- 
ciated with the volume-weighted mixing law is consi- 
dered as a methodology to get in a simple way the com- 
bustion average characteristics. The Eddy-dissipation 


model is used to describe the turbulent combustion. A 
one-step chemical reaction is proposed. The obtained 
results are confronted to the available experimental data. 


Configuration and operating conditions 


The test case RCM-3 presented at the second IWRCM 
workshop [14] is considered as a reference for validation 
(Figure 1). In this test case the single element combustor 
Mascotte V03 is fed with liquid oxygen at 85 K and ga- 
seous hydrogen at 287 K Table 1. The high pressure 
combustion chamber is a square duct with a 50 mm inner 
dimension and with a length of 400 mm. At its down- 
stream end, there is a nozzle of variable shape having a 
convergent length of 20 mm and a throat diameter of 9 
mm allowing a chamber pressure of about 60 bars. 
However, due to the lack of experimental data about the 
nozzle flow, it is replaced by a non-reflecting constant- 
pressure outlet and the cooling film of the windows is 
neglected in the simulation. 
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Fig. 1 View of the Mascotte test-rig (V03) [1] 


Table 1 Operating conditions 


Operating conditions Hə O2 
Pressure [bar] 60 60 
Temperature [K] 287 85 
Mass flow rate [kg.s"] 0.070 0.100 
Velocity inlet [m.s] 236 4.35 
Density at injector inlet [kg/m’] 5.51 1177.8 


Phenomenological trend 


As described in the references [2] and [30], mixing is 
a key point for thrust and efficiency of combustion sys- 
tems. It becomes crucial in the case of liquid rocket en- 
gines because the pressure chamber often exceeds the 
critical point of loaded propellants. Mixing becomes an 
important scientific issue as fluid properties differ from 
classical ideal gas assumption. Liquid and gaseous phas- 
es are no longer separated above the critical pressure and 
the temperature of the mixture. Experimental investiga- 
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tions confirm that surface tension of the oxygen liq- 
uid-core vanishes [42]. The sharp distinction between gas 
and liquid disappears and the entire field becomes essen- 
tially a continuous environment, with no abrupt phase 
change involved in the burning process. On the liquid 
side of the critical isotherm, the material behaves more 
like a very dense gas than a liquid. Near the critical point, 
small state changes have great effects on transport prop- 
erties and variables of state, leading to huge gradients of 
density and other thermodynamic variables during the 
mixing of trans or supercritical fluids [42]. These condi- 
tions make it very difficult to accurately predict the 
processes taking place. Currently some works are effec- 
tuated trying to give new theoretical description that qu- 
antifies the effects of real-fluid thermodynamics on liquid 
fuel injection processes as a function of pressure at typi- 
cal engine operating conditions [21,33,34,38,39,43-45]. 


Governing equations 


In this study, a simplified approach based on stable 
supercritical fluid assumption is adopted. Mass, momen- 
tum, energy and chemical species governing equations 
are written respectively in a conservative form, for a 
compressible and fully turbulent flow. The body forces 
and the radiation effects are neglected. The surface ten- 
sion is also negligible compared to aerodynamic forces 
(infinite Weber number). 


Op = = 
L 4V. =0 1 
wa (1) 
oô( pü P 3 
+ pii(V ii )=-Vp+Vr (2) 
Ot 
where 
2 (2 | = 2 St > 
r=<u(V a)r +(V +V" a) (3) 


t and I are the viscous stress and the identity tensor re- 


spectively and u is the molecular dynamic viscosity; 
O(ph,) oO = 3 Z 
L JOP) a 0n =—pV-u+tc-u-V-g (4) 
t t 


where h, and q are the specific total enthalpy and the 


total heat fluxes, respectively. When neglecting thermal 
diffusion caused by temperature gradient and molecular 
diffusion caused by species concentration gradient (Soret 
& Dufour effects), the heat flux can be expressed as: 


G=-AVT +> pDyhyVYy (5) 
a 
where 

T 
h= ote ae (6) 

T? 
h, is the partial enthalpy of the i" species in the mixture; 

O( py, as 3 r 

Cera), pi WY, =¥-(pD,WY,) +a, (7) 
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©, is the instantaneous rate of production of the i” 


chemical species due to the reaction. 

The Favre decomposition leads to the mean quantities 
of the reactive flow. Every extensive quantity ø (except 
density and pressure), can be decomposed into a mean 
part ø and a fluctuation part g as: 

p=p+o" (8) 

The Favre (mass weighted) averaging form is defined 
as: 

g-— (9) 
p 

Favre decomposition generates constraints due to the 
appearance of fluctuation terms that need to be modeled. 
Hence, the momentum governing equation decomposi- 
tion leads to the following constraint: 

Ox; Ox; 3 Ox, 
for which the selected closed model in this work is the 
(k-€) turbulence model. The other constraints corres- 


ponding to species and enthalpy fluctuations are modeled 
by the following gradient expressions: 


S 


(10) 


n on Ya 
pu, a =A (11) 

Sc, OX; 

— r; by oh 
ush" = — — 12 
pus Pr, Ox, (12) 


where Sc, and Pr, are the turbulent Schmidt and Prandtl 
numbers which are supposed to be identical to treat tur- 
bulent diffusion processes for heat and species in the 
same way. 


Properties evaluation 


Fluid mixture properties at critical regime are difficult 
to obtain and the experimental data on multi-component 
mixtures are particularly scarce. For such regime, proper- 
ties of hydrogen and oxygen differ significantly from 
those at subcritical conditions, so the ideal state equation 
describing the fluid mixture behavior becomes invalid. 
To describe the real gas behavior, many equations of state 
are available such as the SRK [46] or the PR [47] which 
are often used in the supercritical CFD applications. 
More accurate equations of state exist such as BWR [48] 
but they are more complex because they comprise the 
critical compressibility factor as third parameter in addi- 
tion to the critical pressure and critical temperature. R.N. 
Dahms and J.C. Oefelein [45] have conducted a study 
using a modified 32-term BWR equation of state. This 
study provides some new indications on the breakdown 
of the LOx-H, molecular interfaces. The cubic equations 
SRK and PR have some limitations; for example, the 
densities of the fluids are not properly evaluated. This is 
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observable in Figures 2 and 3 where the evolution of 
oxygen and hydrogen densities according to temperature, 
obtained from the NIST library, are confronted to those 
obtained from PR and SRK equation of states. Instead, it 
can be noted in the graphs that the application range of 
the cubic equations is limited to relatively diluted zones. 
In these zones, results obtained with the PR and with the 
SRK equations are superposed on the data retrieved from 
the NIST database [49]. In the dense phase, the results 
are in contrast. PR and SRK equations slightly overesti- 
mate the density of H,. For O, they have an opposite be- 
havior, PR equation overestimates the O, density whereas 
SRK equation underestimates it. To overcome these li- 
mitations, we have estimated the hydrogen and oxygen 
property variations according to temperature at 60 bars 
using a polynomial fitting from the NIST database (see 
Appendix Figures Al to A8) instead of exploiting the 
corrected forms of cubic equations which will be vali- 
dated using NIST data. The description of pure oxygen 
density requires multiple polynomials implemented in 
Ansys-Fluent software as a user defined function (UDF). 
In contrast to subcritical conditions, the rate of combus- 
tion at supercritical conditions is mainly governed by the 
mixing process. Experimental investigations indicate that 
the rate of combustion is governed by the process “mix- 
ing-limited” [50,51]. For the mixing law, the density is 
computed using the Volume Weighted Mixing Law 
(VWML) model as: 

l 


Y, 
du Pi 


The transport properties of the mixture like viscosity, 
thermal conductivity, specific heat and mass diffusion are 
based on a simple mass fraction average of the pure 
component species given by expressions (14), (15), (16) 
and (17) respectively. 


p= (13) 


u=} Yit; (14) 
DI (15) 


cp =} Mens (16) 


J, =-[pDn + Hoy, (17) 
Se, j 
where Dim is the diffusion coefficient for species i in the 
mixture. It is provided from polynomial fitting of hydro- 
gen and oxygen diffusion coefficients (see appendix Fig- 
ures A9 and A10). Figure 4 shows our results concerning 
the density evolution for cold flow using both mixing law 
and cubic equations of state (SRK and PR). It can be no- 
ticed that for both PR and SRK equations of state, result 
of mixture density evolution differ remarkably from the 
VWML in the first portion of axial distance. According to 


the experimental visualizations [51] where the presence 
of the dense core (which supposes a higher density level) 
is clearly illustrated along the axial distance of the com- 
bustion chamber x = 0 to x = 0.1 m, it can be deduced 
that the results of density profile showed by Figure 4 are 
underestimated. However, the VWML provides slightly 
better density results near the injector orifice but remains 
still insufficient after the axial distance x = 0.025 m like 
PR and SRK equations. The asymptotic value of mixture 
density is about 20 kg/m’. 


Combustion model 


The Eddy-Dissipation model is used. Reaction rates 
are assumed to be controlled by the turbulent mixing and 
the Arrhenius chemical kinetic calculations can be 
avoided. Hence, most of the fuels are fast burning. The 
net production rate of species 7 due to reaction r is given 
by the smallest result delivered by equations (18) and 
(19). The chemical reaction rate is governed by the large 
eddy mixing time-scale (k/e). Combustion proceeds 
whenever k/e > 0 and an ignition source is not required to 
initiate combustion. Fluent sets the product mass frac- 
tions to 0.01 which is usually sufficient to start the reac- 
tion. 


E Y 
@,.=v! M_ Ap- min,| ——— 18 
Ir ir @,i PT (| ( ) 
Y 
O; = v! M„;4BpŽ a - (19) 
k È VirM oi 


where Y, is the mass fraction of any product species p, Yr 
is the mass fraction of a particular reactant R, 7’; and 1", , 
are respectively, the stoichiometric coefficient for reactant 
i in the reaction r and the stoichiometric coefficient for 
product 7 in the reaction r. M}; is the molecular weight. A 
and B are empirical constants equal to 4.0 and 0.5 re- 
spectively. 
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Fig. 2 Pure O, density evolution (P = 60 bars) 
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Fig.3 Pure hydrogen density evolution (P = 60 bars) 
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Fig.4 Mixture density (cold flow simulation, P = 60 bars) 
Simulations 


This work is performed with the CFD code Ansys- 
Fluent. The RANS approach is considered due to its sig- 
nificant advantages over the computing time and the pos- 
sibility to simulate large computing domains. The closure 
model (k-e) is selected here due to its advantage of being 
used in harsh conditions of industrial-type flows with 
complex geometries and high Reynolds numbers. Due to 
the chamber symmetry, only one quarter of the whole 
volume is considered in order to reduce the computation- 
al cost. The 3D computing domain is displayed in Figure 5. 


Meshing Method 


The mesh was created so as to have a flagged and ad- 
justed boundary to dispose of the real geometry. The grid 
control and refinement process is achieved using the 
AMR technique in order to improve the quality of the 
results in the zones of strong gradients. The AMR algo- 
rithm creates adaptively finer grids depending on where 
an extra solution is needed in the domain according to a 
user defined error estimator. In the present work, the 
AMR is activated for the mean reaction rate. The original 
cell size (Figure 6a), which is determined before the 
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simulation, starts with cell size of 0.15 mm at the injector 
exit and becomes coarse in the x, y, and z directions with 
the goal of capturing mixing layers (Figure 6b). AMR is 
activated because it is difficult to determine a priori dur- 
ing the simulation where embedding should be added. 
AMR algorithms automatically add the embedding where 
the flow field is the most under-resolved or where the 
sub-grid is the largest. Consequently, the cell size and the 
total number of cells are not constant and change during 
the simulation. For example: the total cell number at the 
start of simulation was 794270 cells, and it became 
915349 cells at the end of the computation. Figure 6c 
shows the final computational grid domain with unstruc- 
tured mesh. 
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Fig.5 Geometrical modeling and boundary conditions 


Boundary Conditions 


In order to solve the previous governing equations, 
boundaries and boundary conditions must be specified 
for each equation. The type of wall boundary is set as 
stationary and smooth boundaries. Law of the wall is 
used for gas temperature and velocity. Both for velocity 
profile of GH, and LOx at the inlet of the combustion 
chamber, preliminary simulations of injectors were 
achieved. This work is done separately with extracting 
the injector to simulate flows in order to deduce the ve- 
locity profiles of GH, and LOx. The inlet and outlet con- 
ditions of the injector flow are imposed respectively ac- 
cording to the injector feed conditions (Table 1) and the 
combustion chamber conditions (pressure). Figure 7 
shows the semi-profile velocity of GH, and LOx along 
the radial direction. The initial turbulence parameters are 
chosen according to those suggested at the second 
IWRCM workshop [14] for the test case RCM-3. 
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Fig. 7 Velocity profile, a. GH, b. Lox 
Results and discussions 


To ensure a good initial computation for the hot flow 
simulation, the first step consists in solving species 
transport equation (7) without considering the chemical 
source term. This resolution is called cold flow simula- 
tion. In this case, the contour of the axial velocity (Figure 
8) obtained using the (k-e) model shows clearly a recir- 
culation zone near the chamber corner due to the pres- 
ence of the wall, and fluid recirculation zone which ex- 
tends to the downstream of the mixing layer. This is con- 
sidered as a consequence of the shearing between the 
propellant flows. 

For hot flow simulation, the combustion model EDM 
predicts the maximum temperature of about 4200 K by 
using a one-step chemical reaction only with H,O as 
combustion product. This value is higher than the hydro- 
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(pore e 


ee 
E nE 


-2.97e-01 


212402 


Fig. 8 Axial velocity field on 45° inclined plane (y = 25 mm, 
z = 25mm) 


gen/oxygen adiabatic flame temperature which is about 
3600 K obtained from the computer program [52]. So, 
the dissociation effect is not taken into account and con- 
sequently the reaction requires additional products. The 
integrating of detailed chemical kinetics model from 
Chemkin [53] significantly adds to the computational 
time. Hence, a two-step chemical reaction is considered 
by adding new species. Several mechanisms have been 
tested [28,53-55]. Among them, the model used by Be- 
narous et al. [35] has provided acceptable results with the 
additional species OH. 
H: + 0.5 O2 > H,O 
H +O, > 20H 

According to this mechanism, the obtained flame 
temperature was improved. Moreover, we have proposed 
and tested the following one-step generic chemical reac- 
tion including hydrogen peroxides (H202) and OH. This 
schema has also provided acceptable results with less 
computational cost: 

2.5H,+ 20, > H20 + OH + H20, 

Figure 9 illustrates the evolution of axial temperature. 
It can be observed that the obtained results using the ge- 
neric reaction are better than those with the two step- 
reactions. It is shown that the axial temperature evolution 
of the flame resulting from EDM model provides a 
maximum temperature at the axial position of ~100 mm 
from the injector and it can be recognized that this posi- 
tion indicating the flame tip is located at about the same 
axial location (~100 mm) according to the temperature 
measurements [51]. 

The experimental data [1] indicate that typical 
GH,/LOx flames are directly attached at the tip of the 
LOx post in contrast to what is generally observed for 
hydrocarbon flames which are lifted-off and the general 
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structure of the flame can be characterized as a thin 
sheath around the dense gaseous oxygen. This intensive 
reaction zone results from the absence of atomization for 
a supercritical fluid, which leads to the occurrence of 
only thin mixing layers between LOx and GH). As a 
consequence, even for fully turbulent flow cases, the 
flame thickness not far from the orifice remains limited 
and its front is confined close to the LOx jet. The flames 
resulting from numerical simulations using the generic 
one-step chemical reaction and the two-step reaction 
model are shown in Figure 10. It is clear that the flame is 
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Fig.9 Axial temperature evolution 
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(a) Generic one-step chemical reaction with standard (k-e) 
turbulence model 


(b) Generic one-step chemical reaction with realizable (k-s) 
turbulence model 


(c) Two-step reaction with standard (k-s) 
turbulence model 


(d) Two-step reaction model with realizable (k-e) 
turbulence model 
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Fig. 10 Contour of temperature 
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stretched with the use of the (k-e) realizable turbulence 
model. The flame obtained using the standard (k-£) tur- 
bulence model is better and with a realistic flame length. 
Figure 11 confronts the obtained OH species distribution 
with the experimental visualization. 

Along the axial direction, the density and the dense core 
were differently reported in experimental investigations. 
Flow visualization indicates that at high pressure, the 
LOx jet is no more purely liquid, but becomes a super- 
critical fluid, even if it remains very cold and dense [56], 
but more accurate information about the axial density is 
still difficult to capture. Figure 12 depicts axial density as 
a function of temperature and Figure 13 displays the axi- 
al evolution of the mixture density both for cold and hot 
flow. The cold density profile shows a two-step density 
evolution, rapid and slow decrease, indicating the transi- 
tion of oxygen from transcritical to supercritical regime 
close to the injector tip. The hot flow density shows a 
progressive decreasing curve, which could be explained 
by the presence of the flame sheath which prevents cro- 
ssing oxygen gas or ligament. It can be observed that the 
hot flow density seems to be incorrect far from the injec- 
tor which indicates the limit of the VWML approach. 

Figure 14 shows a qualitative comparison of the radial 
density with available experimental investigations [51] 
using the light intensity transmitted by the flow at 
x/Drox=0.2. The light intensity transmitted by the LO x jet 
can exhibit density gradients. For dark location, the light 
absorption is strong and the intensity remains small indi- 
cating the presence of the dense core. For gaseous phase, 
intensity becomes higher for small densities. 


Fig. 11 OH concentration. (a) Experimental visualization [37]. 
(b) Present simulation (generic one-step chemical 
reaction with standard (k-s) model) 
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Fig. 12 Mixture density confronted to pure hydrogen and 
oxygen density along the centerline 
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Fig.13 Mixture density comparison 
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Fig. 14 Density at the position x/D;,,, = 0.2 confronted to 
measured light intensity 


Conclusions 


A 3D RANS study of LO 2/GH) supercritical combus- 
tion is performed using the one-phase multi-component 
approach on the Mascotte test-rig (V03) case. Polynomial 
fitting of the pure components O, and H, properties has 
been accurately accomplished with the NIST library da- 
tabase, especially close to the critical temperature since 
propellant properties have an important influence on re- 
sults. These polynomials are implemented in Ansys- 
Fluent software as a user defined function. The numerical 
study is carried out in two steps. Firstly, the cold flow 
simulation is realized to appreciate the quality of the tur- 
bulence model and to get best starting conditions for hot 
flow computations. The hot flow simulation is achieved 
using the EDM combustion model associated with two 
chemical reaction models. The first one is the two-step 
chemical reaction scheme already used by [35] in a 2D 
simulation and the second one proposed in this work is a 
generic one-step chemical reaction. Compared to the cu- 
bic equations of state, the obtained results indicate a bet- 
ter estimation of the mixture density close to the injection 
orifice when the volume-weighted mixing law is used. 
Results concerning centerline temperature are in good 


agreement with experimental data and the OH concentra- 
tion has a good qualitative concordance with the ob- 
served one. The radial evolution of the density at the non- 
dimensional axial station x/D;,,= 0.2 where experimental 
data are available also seems to have an acceptable evo- 
lution. However, much work needs to be carried out in 
order to better capture the hot mixture density in the 
whole combustion chamber by improving the mixture 
rule (VWM.L) and/or the cubic equation of state. 
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Appendix A 
Properties of pure hydrogen and oxygen presented at the pressure of 60 bars. Figures Al to A8 show density, specific 
heat, thermal conductivity and viscosity evolutions (NIST Chemistry WebBook). Figures A9 and A10 show the diffu- 


sivity evolution [21]. 
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Fig. Al Hydrogen density at constant pressure [60 bars] Fig. A2 Oxygen density at constant pressure [60 bars] 
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Fig. A3 Hydrogen specific heat at constant pressure [60 bars] 
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Fig. A4 Oxygen specific heat at constant pressure [60 bars] 
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Fig. A5 Hydrogen thermal conductivity at constant pressure [60 bars] 
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Fig. A6 Oxygen thermal conductivity at constant pressure [60 bars] 
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Fig. A7 Hydrogen viscosity at constant pressure [60 bars] Fig. A8 Oxygen viscosity at constant pressure [60 bars] 
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Fig. A9 Hydrogen diffusivity at constant pressure [60 bars] 
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Fig. A10 Oxygen diffusivity at constant pressure [60 bars] 


